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§1 .  Introduction 

In  a  previous  paper  (Holmes  and  Marsden  [1981b])  we  developed  a  method 
for  proving  the  existence  of  Smale  horseshoes  in  two  degree  of  freedom 
Hamiltonian  and  nearly  Hamiltonian  systems.  This  paper  extends  those  methods 
to  systems  with  three  or  more  degrees  of  freedom.  We  start  with  an  unper¬ 
turbed  system  containing  a  homoclinic  orbit  and  at  least  two  families  of 
periodic  orbits  associated  with  action  angle  coordinates.  We  use  KAM  theory 
to  show  that  some  of  the  resulting  tori  persist  under  small  perturbations 
and  use  a  vector  of  Melnikov  integrals  to  show  that,  under  suitable  hypotheses, 
their  stable  and  unstable  manifolds  intersect  transversally.  This  transverse 
intersection  is  ultimately  responsible  for  Arnold  diffusion  on  each  energy 
surface. 

Our  methods  are  a  generalization  of  those  of  Arnold  [1964]  where  "Arnold 
diffusion"  was  first  introduced.  The  applications  are,  however,  somewhat 
different  and  we  believe,  of  more  direct  physical  interest. 

The  main  example  treated  in  this  paper  is  to  a  Hamiltonian  system  con¬ 
sisting  of  a  pendulum  coupled  to  two  oscillators  (with  amplitude-dependent 
frequencies).  The  system  is  shown  to  have  Arnold  diffusion.  Using  the 
techniques  in  our  previous  paper,  one  can  also  show  that  the  Arnold  diffusion 
on  a  certain  energy  surface  survives  suitable  positive  and  negative  damping 
perturbations. 

We  shall  assume  that  our  coordinates  are  given  in  canonical  form.  How¬ 
ever,  many  interesting  problems  involving  rigid  body  dynamics  are  best  done 
in  a  more  general  Lie  group  theoretic  context.  This  situation  is  discussed 
in  Holmes  and  Marsden  [1981c]. 

Acknowledgements.  We  thank  Alan  Weinstein  for  several  helpful  discussions  and 
A1  lan  Kaufman  for  suggesting  a  stimulating  physical  example. 
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§2.  Transversal  Intersection  of  Invariant  Manifolds  by  Melnikov's  Method 


In  this  section  we  are  concerned  with  perturbations  of  Hamiltonian 
systems  of  the  form 


H  (q.p.x.y)  =  F(q,p)  +  G(x,y)  » 


(2.1) 


where  (q,p,x,y)  are  canonical  coordinates  on  a  2(n  +  1)  dimensional 
symplectic  manifold  P;  q  and  p  are  real  and  x  ~(x\  ....  xn), 
y  =  (y1 ,  ....  yn).  We  assume  the  coordinates  are  canonical,  although,  in 
some  examples  such  as  the  rigid  body  this  requires  modification  (Holmes  and 
Marsden  [1981c]).  We  shall  also  assume  that  action  angle  coordinates 
(Q-j ,  •••,  9n,  1-j,  ...»  In)  can  be  found  in  a  certain  region  of  phase  space 
such  that  (2.1)  takes  the  form 


H°(q ,p ,9 ,1 )  «  F(q,p)  +  £  G.(I.)  . 

i=l  1  1 

We  also  assume  that  G(0,0,  . . . ,  0)  *  0  and  that 


(2.2) 


W-3it>0  ,or  V° 

J 


(2.3) 


The  perturbed  problem  we  consider  has  the  form 


n  1 

He(q,p,0,I)  =  F(q.p)  +  £  G  (I)  +  eH*(q.P.e,I) 

i=l 


(2.4) 


where  H  is  2^-periodic  in  9^,  ...,  9n.  Now  we  recall  how  this  n+1  degree 
of  freedom  system  may  be  reduced  to  an  n  degree  of  freedom  non-autonomous 
system;  the  reader  should  refer  to  Holmes  and  Marsden  [1981b]  for  details. 

Choose  one  of  the  action  coordinates,  say  In-  Since  ('!,,)  >  0  for 
In  >  0,  we  can  invert  the  equation 


*— t  -v«.  - 


fit*'' 


I 
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a 


i1 

ti 


H£(q,p,e,I)  =  h. 


to  obtain 


If  we  write 


-  L  (q,P.  0-] »  •••»  ©n»  1 1 »  ...»  In_i ;  h)  . 


L£  -  L°  +  eL1  +  0(e2)  , 


(2.5) 


(2.6) 


(2.7) 


then  a  simple  computation  shows  that 


L  (q,P,  1 1 »  ...»  I^_-j  >  h) 

.1  n-1 

=  G  ‘(h  -  F(q,p)  -  l  G  (I  ))  , 

j=l  J  J 


(2.8) 


L  Bn »  I  -j  »  •••»  I  ^  1>  h ) 

H  (q,p,0^»  ...,  0n»  I-j  »  ...»  In_i  ,  L°(qj,p,I-|.  ;h) 

^n(L  (q,P,Ii »  ...»  I  i»  h) 

(2.9) 

Changing  variables  from  t  to  9n  and  writing  (  )'  for  (  ),  Hamilton's 

,  c  n 

equations  for  H  become 


s  3Lf.  n.  .  ikl 

3p  *  P  ~  3q 


A,  _  5Ll  T,  _  3L£  .  . 

j  '  ai,  ’  j  ”  ‘  ae,’  J  "  *•••  n_1 

J  J 


(2.10) 


Using  (2.7)  -  (2.9),  equations  (2.10)  are  in  the  form  of  a  2iT-periodically 
forced  n  degree  of  freedom  Hamiltonian  system.  Notice  that  L°  is  (formally) 
completely  integrable,  having  n  constants  of  the  motion  given  by 


or,  al ternatively. 


L°(energy)  and(l],  ....  1^)  *  (J^,  ....  H^) 

L  and  (G-j(I^)»  ...»  Gn_^(I^_^))  *  (h^,  ...,  hn_.|) 


4 


(This  reflects  the  general  fact  that  complete  integrabil ity  is  preserved  by 
the  reduction  process). 


Assume  now  that  the  Hamiltonian  F  has  a  homoclinic  orbit  (q(t),  fT(t)) 
joining  a  saddle  point  (qo,pQ)  to  itself.  [The  case  of  heteroclinic  orbits 
connecting  different  saddle  points  proceeds  in  the  same  way].  The  Hamiltonian 
system  for  L°  thus  has  an  n-1  parameter  family  of  invariant  n-1  dimensional 
tori  T(h^,  ...,  hn  given  by 

G.(I.)  =  h.  =  constant  (i.e.  I.  =  l.  -  G.^h  )) 

J  J  J  J  W  J  J 

0j  -  +  9j(°)  Cmod  2tt),  j  =  1,  ....  n-1  (2.11) 

q  *  q0»  P  =  P0 


Correspondingly,  the  system  for  F  has  an  n  parameter  family  of  invariant 
tori  T(h1>  ....  h^).  Henceforth  we  write  the  (phase)  constants  of  integra¬ 
tion  9.(0)  as  9°,  j  =  1 ,  . . . ,  n-1 ,  n. 

J  J 

The  torus  T(h^,  ...»  h  .j )  is  connected  to  itself  by  the  n-dimensional 
homoclinic  manifold 


6j  J  =  ' . • 

q  *  q(9  -  9°).  P  *  p(6  -  9°)  , 

where  the  phase  constant  9°  associated  with  the  'reduced*  degree  of  freedom 
appears  explicitly.  This  manifold  consists  of  the  coincident  stable  and  un¬ 
stable  manifolds  of  the  torus  T(h^,  ...,  hn_-j ) ;  i.e. 


Ws(T(hr 

given  by  (2.12).  See  Figure  1. 


V,))  -  Wu(T(hr  ...,  hf|_1 ) ) 


•  •  •  * 
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I .  =  hi ,  0 .  27r-periodic, 
J  J  J 

X 


T (h-j »  •  •  •  >  ) 

Figure  1 

For  e  f  0  the  system  (2.10)  possesses  a  Poincare  map  from  (a  piece 

of)  (q,p,  0-j ,  ....  9n_i»  !-|  >  •••»  In_i)  space  to  itself  where  0n  goes 

through  an  increment  of  2n,  starting  at  some  fixed  value  0°,  (which  will  be 

suppressed  in  the  notation).  The  tori  T(h.,  ....  h  .)  are  invariant  mani- 

i  n- 1 

folds  for  P  .  In  fact,  these  tori  are  isotropic  submanifolds  (i.e.  the 
canonical  2-form  u  vanishes  on  them),  a  fact  we  shall  need  later. 

The  program  is  to  show  that  for  e  f  0  some  of  the  tori  persist  and 
that  their  stable  and  untable  manifolds  intersect  transversely.  To  do  this 
we  shall  invoke  the  KAM  theory  and  Melnikov's  method.  The  result  will  then 
be  interpreted  as  Arnold  diffusion. 

Let  us  first  discuss  the  invariant  tori.  The  manifold  obtained  by  setting 
q  =  q0,  P  =  PQ  is  a  2n-2  dimensional  normally  hyperbolic  invariant  manifold, 
say  Mq ,  for  our  Poincare  map  P  .  Thus,  for  e  small,  MQ  perturbs  uni¬ 
quely  to  an  invariant  manifold  M£  for  P£.  The  KAM  thoery  now  can  be  applied 
to  the  family  of  invariant  tori  T(h^,  ...,  hn  ^)  on  M  .  If  the  hypotheses 
of  non-degeneracy  and  non-resonance  hold  (see  Arnold  [1978,  Appendix  8])  then 

the  torus  T(h,,  ....  h  .)  will  perturb  to  an  invariant  torus  T  (h.,  .... 

•  n- 1  el 

h  j )  for  P  ,  for  e  sufficiently  small  (depending  upon  the  precise  'degree' 
of  nondegeneracy).  Moreover,  the  proof  shows  that  this  torus  is  also  isotropic 
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(see  Moser  [1973]).  We  note  that  the  perturbed  torus  T£  has  the  same  fre¬ 
quencies  n.(h.)  as  the  unperturbed  torus  and  thus  the  perturbed  phase  angles 

J  J 

do  not  drift  appreciably  from  the  unperturbed  ones.  We  use  this  fact  below. 

Although  a  set  of  positive  measure  of  the  perturbed  tori  persist  near 
the  original  ones,  the  resonant  tori  containing  continuous  families  of 
periodic  motions  generally  break  into  finite  sets  of  alternating  elliptic 
and  hyperbolic  periodic  orbits  with  associated  homoclinic  motions,  as  in 
Arnold  [1978],  p.  397.  The  boundaries  of  the  elliptic  islands  are  conven¬ 
tionally  drawn  as  homoclinic  orbits  of  a  flow:  these  actually  belong  to  an 
associated  averaged  (=  canonically  transformed)  system.  Restoration  of  the 
terms  omitted  in  averaging  leads  to  the  prediction  that  these  islands  will, 
in  turn,  be  surrounded  by  regions  containing  transverse  homoclinic  orbits, 

(cf.  Holmes  [1980])  but  these  regions  are  smaller  than  any  power  of  e, 
since  they  can  be  removed  by  successive  averaging  operations.  In  fact  such 
'stochastic  layers'  are  generally  exponentially  small  in  e  and  attempts  to 
compute  them  by  the  Melnikov  method  necessitate  a  careful  examination  of 
errors.  This  will  be  the  subject  of  a  further  publication;  cf.  Sanders  [1980]. 

In  the  case  of  two  degrees  of  freedom  for  which  the  unperturbed  reduced 
system  has  a  hyperbolic  saddle  point  xQ  =  (90*P0)»  solutions  of  the  per¬ 
turbed  system  lying  in  the  perturbed  stable  and  unstable  manifolds  of  the 
perturbed  saddle  point  x£  of  the  map  P£  can  be  expanded  in  power  series 
which  converge  uniformly  in  the  intervals  indicated: 

WS  :x*  =  x(9  -  9°)  +  ex?(8,0°)  +  0(e2),  0  £  [0°,»), 

e  1  (2.13) 

WU  :  x^  =  x(9  -  9°)  +  ex^e.e0)  +  Cl(c2),  0  €  (-»,9°], 
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%  I 

iV 

»  w 

y 


o  h(e  ‘  e°M 

where  x(9  -  0  )  =  is  the  unperturbed  homoclinic  orbit.  (Recall 

P(e  -  e°) 

that  the  periodic  variable  0  has  replaced  time.)  For  details,  see  Holmes 
[1980],  Sanders  [1980]  or  Greenspan  and  Holmes  [1981].  [Basically  (2.13) 
follows  from  the  fact  that  the  perturbed  solutions  lie  in  manifolds  of  solu¬ 
tions  forward  or  backward  asymptotic  to  the  perturbed  saddle  points.  )  Similarly, 
solutions  lying  in  the  perturbed  invariant  manifolds  WS(T£),  wu(Te)  of  a  per¬ 
turbed  torus  T£  can  be  expanded  in  convergent  power  series  in  e  in  such 
intervals,  since  the  perturbed  actions  are  e  close  and  the  perturbed  angles 
do  not  drift  but  remain  close  to  the  unperturbed  angles  on  the  tori.  This 
result  will  be  used  implicitly  in  what  follows. 

The  perturbed  invariant  manifolds  WS(T£)  and  WU(T£)  of  the  torus 
T  for  the  map  P  are  n-dimensional  manifolds  lying  Cr  close  to  the  un- 
perturbed  homoclinic  manifold  given  by  (2.12).  i.e. 


F  s  h,  1.  s  l,  j  5  1,  ....  n-1 

J  J 


(2.14) 


where  h  is  the  energy  of  the  homoclinic  orbit  for  F.  Now  we  are  ready  to 
give  a  criterion  for  the  transversal  intersection  of  WS(T£)  and  WU(T£).  In 
order  for  the  results  to  be  applicable,  it  is  useful  to  present  the  hypotheses 
in  terms  of  data  given  for  the  original,  rather  than  the  reduced  system. 

We  consider  a  Hamiltonian  system  with  n+1  (>^  3)  degrees  of  freedom 
of  the  form 


He(q,p,  0lt  ...,  6n,  Ir  ....  IR)  =  F(q,p)  +  J  G^L) 


H  (q,p,0.j*  •  •  • »  9  ^  j  1 1  >  •••*  I  ^ ) » 


(2.15) 


Introduce  the  following  assumptions  and  terminology 
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(HI)  F  contains  a  homoclinic  orbit  (q(t),  p(t))  connecting 
a  saddle  point  (q0,pQ)  to  itself.  Let  F  be  the  energy 
of  this  orbit. 

(H2)  £2j(Ij)  =  Gj ( I j )  >0  for  j  -  1 . n. 


^4 


4 


L 


Let  h  >  h  and  let  the  unperturbed  homoclinic  manifold  be  filled  with 
an  n-parameter  family  of  orbits  given  by  (F>P>  ^ »  • ..,  6n,  1^ ,  ...,  I  )  = 
(q(t),p(t),  Sij ( 1^ )t  +  9°,  ...,  ^n(In)t  +  9°,  I-j ,  ....  In).  Pick  one  such 
orbit  and  let  {F,H^ }  denote  the  (qtp)  Poisson  bracket  of  F(q,p)  and 

H\q,p,  0^,  ....  0n>  1^,  _ ,  I  )  evaluated  on  this  orbit.  Similarly,  let 

{Ik>H^}  =  -3HV30k»  k  =  1,  ...,  n-1  be  evaluated  on  this  orbit.  Define 
the  Melnikov  Vector  M (6°)  =  (M^ ,  . . . ,  ^  ,  M  )  by 


Hk(0“, 


Mn<6?' 


9n>  h,  h^ ,  h^, 


V  h,  hr  h2. 


nn  -i )  =  (I.  ,  H  )  dt,  k  =  1 ,  . . . ,  n-1 


Vi}  - 


{F,H  }  dt. 


(2.16) 


(We  note  that  h  =  h  -  h  -  Y  h.;  I  and  h  do  not  explicitly  enter 
n  ^“l  j  n  n 

the  calcualtions,  since  In  is  eliminated  by  the  reduction  process;  we  also 
note  that  these  integrals  need  not  be  absolutely  convergent,  but  we  do  require 
conditional  convergence.) 

(H3)  Assume  that  the  constants  G.(I.)  =  h.,  j  =  1,  ....  n  are 

J  J  J 

chosen  so  that  the  unperturbed  frequencies  £2,(1,),  ...,  £2,(1  )  satisfy 

II  In 

the  non-degeneracy  conditions  (i.e  £2'.(I.)  t  0,  j  *  1 ,  . . . ,  n-1)  and  the 

J  J 

non-resonance  conditions  mentioned  above  (cf.  Arnold  [1978],  Appendix  8). 


•  v « ' 


i  i , 


(H4)  Assume  that  the  multiply  2tt  periodic  Melnikov  vector  M  :  IRn  IRn 
(which  is  independent  of  e)  has  at  least  one  transversal  zero;  i.e.  there 
is  a  point  (8°,  ....  9°)  for  which 

M(8°,  ....  0°)  -  0 

bUt  det[DM(8|,  ....  0°)]  f  0  , 

where  DM  is  the  n  *  n  matrix  of  partial  derivatives  of  M-j,  ...»  M  with 

respect  to  0°,  - -  6^,  the  initial  phases  of  the  orbit. 

Here  is  our  main  theoretical  result. 


2.1  Theorem.  If  conditions  (HI)  -(H4)  hold  for  the  system  (2.15),  then,  for 
e  sufficiently  small,  theperturbed  stable  and  unstable  manifolds  WS(T£), 
and  Wu(T£)  of  the  perturbed  torus  T£  intersect  transversely.  (See  Figure  2.) 


WS(T  )  nwu(T  ) 
e  e 


Figure  2.  The  stable  and  unstable  manifolds  of  the  invariant  torus  Te 
for  the  Poincar^  map  in  the  reduced  space  for  a  system  with 
three  degrees  of  freedom.  This  figure  occurs  in  (q,p,6,  ,1.) 
space  (one  dimension  (Ii)  is  suppressed),  in  a 
8j  8  fixed  cross-section  for  fixed  total  energy  h. 


tva* 
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Remark.  The  conclusions  imply  that  the  perturbed  system  has  no  analytic 
integrals  other  than  the  total  energy  H6-  and,  for  n  +  1  ^  3,  that  Arnold 
diffusion  occurs.  This  is  discussed  in  the  next  section. 


Proof.  First  we  notice  that  the  brackets  of  the  original  functions  project 
to  corresponding  brackets  of  the  reduced  system: 


aV1)  = 


(F,h'} 


{Ik’L  } 


n^rrr ^k* H  **  k  =  1 »  •  •  •  • n  ■ 1 


(2.17) 


(2.18) 


(see  Holmes  and  Marsden  [1981b],  Prop.  3.1). 

We  next  wish  to  relate  these  brackets  to  a  vector  measuring  the  distance 

between  the  perturbed  stable  and  unstable  manifolds. 

Consider  the  suspended  system  in  (q,p,9-|,  ...,  0n_-j »  I^»  ...»  In  ^ ) 

space.  Pick  a  transversal,  l  ,  to  the  unperturbed  homoclinic  manifold 
s  u  0 

W  (T)  =  W  (T)  in  (q,p,  9-j ,  ....  9  1^ ,  ...,  In  ^ )  space  at  the  point 

(q(0),  p(0),  0°,  ....  6°_1 ,  1^,  ...,  and  at  "time"  6°.  Now  for  t 

sufficiently  small,  WS(T  )  and  WU(T  )  intersect  E  in  unique  points 

e  e  0o 

in  (q,p,  01 ,  ...,  Qn_i ,  1^,  ...,  I  i)  space,  which  we  denote 

x\(e°,  0°)  and  x“(9°,  9°) 

The  unique  trajectories  in  (q,p,  91 ,  - -  0^,  Ij,  - -  Ifj_1  ,  en)  space 

with  these  points  as  initial  conditions  and  "time"  9n  will  be  denoted 

Xs£(e°,  9n)  and  x“(9°,  e„)  . 


"If+f 


■t 
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As  in  Holmes  and  Marsden  [1981a],  a  measure  of  the  distance  between  these 
vectors  and  the  tangent  to  WS(T)  =  WU(T)  in  the  "0n  direction",  i.e. 
the  direction  of  X  Q ,  the  Hamiltonian  vectorfield  of  the  unperturbed  dyna¬ 
mics,  is  provided  by  the  symplectic  form.  Let 


A  (6  ,6°)  =  w(X 
e,n  n  n'  , 


xs  .  vUl  —  A+ 
o’  e  xe'  e,n 


A"e,n  +  0(e  }  • 


Now  as  6n  -*■  +®,  x*  -*•  T£  and  as  0n  -*•  x^  T  ,  so  as  in  Holmes  and 
Marsden  [1981a]  Lemma  5,  we  obtain 


A  (0°)  =  A  (8°,0°)  = 

e,nv  n'  £,nv  n’  n 


eM  (0.°, 
n  1 


n-1  * 


[«(!„)]' 


e°)  2 

—  +  0(e2}. 


(2. 


Note  that  the  integrals  (2.16)  are  well  defined  since  one  integrates  forward 
along  tne  stable  manifold  and  backward  along  the  unstable  manifold  (cf.  eq. 
(2.13)). 

A  crucial  feature  of  this  calculation  is  the  fact  that  A*  n(0n>0°)  -*■  0 

as  0n  -*■  ±°°  since  ui  vanishes  identically  on  T.  This  holds  as  follows. 

The  invariant  tori  are  isotropic  and  A+  (0  ,0°)  *  ui(X  ex?)  where  X  „ 

e,n  n  n'  L°  L° 

and  x*  are  evaluated  on  the  unperturbed  homoclinic  manifold  and  x^  is 

the  solution  of  the  first  variation  equation.  Since  X  is  tangent  to  T, 

5  I- 

x^  necessarily  approaches  a  tangent  to  T,  so  as  T  is  isotropic, 

A+  ■+  0.  We  note  that  in  this  context  the  perturbed  torus  T  may  move 
e,n  r  e 

(by  0(e))  and  need  not  remain  fixed  as  in  the  special  case  treated  by  Arnold 
[1964]  or  as  in  Melnikov's  [1963]  paper  (cf.  Holmes  [1980]). 

Thus,  Mn(9°»  - -  »  0°)  measures  the  leading  nontrivial  component 

of  the  distance  between  WS(T£)  and  WU(T£)  (up  to  a  constant)  in  a  direc¬ 
tion  transverse  to  the  "dynamic"  variable  0(|.  Likewise,  M^(i  *  1,  ..., 
n-1)  measures  the  distance  between  WS(T£)  and  WU(T£)  in  the  direction 


n 
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transverse  to  the  generator  of  the  9^  variable.  The  theorem  now  follows 
from  these  facts.  ■ 

Remarks  1.  One  can  of  course  permute  which  of  the  action  angle  variables 
is  used  for  the  reduction  procedure.  The  remaining  oscillators  must  satisfy 
the  KAM  non-resonance  and  non-degeneracy  conditions. 

2.  J.  Gruendler  [1981]  has  treated  the  2n  dimensional,  periodi¬ 
cally  forced  case  in  which  one  also  has  an  n-parameter  family  of  unperturbed 
homoclinic  orbits  but  in  which  they  are  homoclinic  orbits  to  a  hyperbolic 
saddle  point  x  and  dim  Ws(x)  =  dim  Wu(x)  =  n.  Again  one  obtains  a  gen¬ 
eralized  n-vector  of  Melnikov  functions  each  depending  upon  n  arguments, 
one  of  which  is  the  section  'time'  (0°)  and  the  remaining  n-1  of  which 
serve  to  parametrize  the  family  of  orbits.  The  manifolds  Ws(x)  and  Wu(x) 
are  both  necessarily  isotropic,  so  one  can  proceed  in  a  way  analogous  to 
that  here.  However  no  KAM  theory  is  needed  and  ordinary  horseshoes  rather 
than  Arnold  diffusion  occur.  Gruendler  applies  applies  the  theory  to  the 
case  of  a  periodically  forced  spherical  pendulum. 

3.  The  theorem  can  be  somewhat  generalized.  For  example,  many  inte  - 

grable systems  do  not  decompose  precisely  as  assumed  in  the  form  F(q,p)  + 
n 

l  G.(I.)  and  one  sometimes  finds  that  the  unperturbed  'frequencies1, 
i=i  J  J 


also  depend  upon  (q,p).  If  this  occurs,  and  ft.(I.,q,p)  is  not  constant 

J  J 

on  the  unperturbed  manifold,  then  it  must  be  incorporated  into  the  Poisson 
brackets  (cf.  eq.  (2.17)-(2.13)).  This  situation  will  be  dealt  with  in 
Holmes  and  Marsden  [1981c]. 
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4.  Alan  Weinstein  has  pointed  out  that  even  without  hypothesis  (H4), 


the  stable  and  unstable  manifolds  of  the  perturbed  torus  T£  must  inter¬ 


sect.  This  comes  about  as  follows.  As  in  the  standard  Melnikov  Analysis 


(Holmes  and  Marsden  [1981a]),  pick  a  2n-dimensional  cross  section  Z  .  The 


stable  and  unstable  manifolds  of  T£  for  the  associated  Poincarg  map. 


WS(T  )  and  WU(T  )  are  Lagrangian  submanifolds  of  Z  which  are  coincident 


6" 


at  e  =  0.  Lagrangian  intersection  theory  (Arnold  [1965]  and  Weinstein  [1973]) 
shows  that  the  perturbed  manifolds  must  intersect.  This  observation  gener¬ 
alizes  one  of  McGehee  and  Meyer  [1974].  It  follows  that  the  hypothesis  (H4) 
holds  for  generic  perturbation  terms  .  However  condition  (H4)  allows  one 
to  check  transversal ity  in  specific  cases. 

5.  In  contrast  to  our  results,  Easton  and  McGehee  [1979]  use 
Moser's  [1978]  fixed  point  theorem  to  show  that  some  homoclinic  orbits  in  a  r>\cdti  Syjte: 
survive  under  special  perturbations.  Alan  Weinstein  points  out  that,  similarly  ? 


at  least  two  homoclinic  orbits  survive  perturbations  of  the  spherical 


pendulum's  S1  family  of  homoclinic  orbits. 


1 
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§3.  Non-integrabil  ity  and  Arnold  Diffusion 

If  the  stable  and  unstable  manifolds  WS(A),  WU(A)  of  a  hyperbolic 

invariant  set  A  intersect  transversely  then  it  follows  fronr  the  lambda 

lernna  (Palis  [1969],  Newhouse  [1980])  that  WS(A)  accumulates  on  itself 

and  WU(A)  accumulates  on  itself.  A  similar  result  holds  for  the  invariant 

tori  in  the  present  case  (cf.  Arnold  [1964],  Theorem  1);  more  precisely,  if 

A-|  c  Wu(T£)  is  an  n  dimensional  neighborhood  of  a  transverse  homoclinic 

point  x  G  (T)?iws(Te) ,  and  C  WU(T£)  is  any  open  disc  then  there 

are  points  of  U  p”(A,  )  lying  arbitrarily  close  to  A0.  Such  a 
n*o  E  1  c 

torus  T£  is  said  to  be  a  transition  torus.  The  torus  is  said  to  lie  in 

.  12  i ( 

a  transition  chain  of  transition  tori  T_,  T  ,  ...  T  if  the  unstable  mani- 

-  £  £  £ 

fold  WU(T^)  of  the  jth  torus  transversely  intersects  the  stable  manifold 

of  the  (j  +  l)st.  This  holds  in  our  case,  since,  by  KAM  theory,  the  set 

of  'sufficiently  irrational'  tori  preserved  when  e  f  0  has  measure 

u(e)  -*-1  as  e  -*■  0  (it  is,  in  fact,  a  Cantor  set).  Thus,  for  sufficiently 

small  e  one  can  find  tori  T^,  T^  which  are,  along  with  their  stable 

and  unstable  manifolds,  arbitrarily  Cr  close  away  from  the  torus  int! 

wrifeb have  large  "oscillations"  near  the  torus  as  in  Holmes  and  Marsden 

[1981b],  Fig.  B.l.  It  follows  that  if  WU(T^)  intersects  WS(T^)  trans- 

£  £ 

versely  it  must  also  intersect  WS(T^)  transversely.  Applying  the  same 
argument  to  T^  ,  T^,  ...  one  constructs  a  transition  chain.  Orbits 
lying  in  WU(TJ)  therefore  accumulate  on  Wu(Tk)  for  k  >_ j  and  these 
orbits  and  nearby  ones  provide  a  mechanism  by  which  solutions  can  "diffuse" 
from  the  neighborhood  of  the  torus  to  any  other  in  the  transition  chain. 

(cf.  Arnold  [1964],  Theorem  2).  An  argument  analogous  to  that  above  shows 
that  WU(T^+1 )  intersects  WS(T^)  and  thus  that  diffusion  can  take  place 


in  both  directions  along  the  chain.  Note,  however,  that  the  length  of  the 
chain  is  generally  governed  by  the  perturbation  strength  e,  since  as  e 
increases  the  set  of  perturbed  tori  generally  diminishes. 

The  mechanism  outlined  above,  which  we  attempt  to  portray  in  Figure  3, 
is  the  basis  for  Arnold  diffusion.  Clearly  it  can  only  occur  in  systems 
with  three  or  more  degrees  of  freedom  (n  >_  2),  since  the  unperturbed  2n 
dimensional  reduced  Poincard  map  must  admit  continuous  families  of  tori 
connected  by  smooth  homoclinic  manifolds, and  this  cannot  occur  in  two 
dimensions.  For  more  information,  numerical  examples  and  physical  insights, 
see  Chirikov  [1979]  and  Lieberman  1.1980].  The  main  physical  consequence  of 
diffusion  is  that  (given  sufficient  time)  energy  can  be  transferred  back 
and  forth  in  relatively  large  amounts  between  distinct  physical  components 
or  vibration  modes  of  the  system.  Moreover  this  transfer  of  energy  will 
typically  take  place  in  an  irregular  manner,  in  contrast  to  the  regular 
quasi  periodic  energy  transfer  occuring  between  modes  in  linear  or  other 
integrable  systems. 

Thus,  in  contrast  to  the  two  degree  of  freedom  case,  in  which  the 
sufficiently  irrational  invariant  tori,  preserved  for  small  perturbations, 
serve  as  boundaries  to  regions  of  homoclinic  (chaotic) motions  in  the  three 
dimensional  total  energy  manifold,  in  systems  with  three  or  more  degrees  of 
freedom  the  solutions  can  diffuse  from  torus  to  torus  along  transition  chains 
the  n-tori  which  are  preserved  do  not  bound  regions  of  2n  +  1  space  for 
n  _>  2.  In  our  case,  since  two  way  transition  chains  can  be  chosen,  we  can 
find  periodic  motions  of  arbitrarily  high  period  close  to  such  chains,  just 
as  in  the  standard  two  dimensional  horseshoe  example.  The  density  of  the 
set  of  such  motions  and  the  dense  orbit  accompanying  them  guarantees  non- 
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Figure  3.  Intersections  of  manifolds  and  Arnold  diffusion 

in  a  three  degree  of  freedom  system.  The  Poincarg 
section  82  =  is  sil0wn  on  the  energy  surface  He  =  h. 


existence  of  any  additional  analytic  integrals  other  than  the  total  energy 
He.  (In  fact  one  sees  that  such  dense  orbits  exist  within  neighborhoods  of  any 
transverse  homoclinic  orbits  connecting  a  torus  to  itself,  without  invoking 
the  idea  of  diffusion.) 

The  presence  of  a  small  amount  of  noise  in  a  system  is  believed  to 
"stabilize"  in  some  sense  the  occurrence  of  Arnold  diffusion,  in  the  same 


way  that  noise  often  "stabilizes"  or  "makes  visible"  horseshoes  (cf.  Holmes 
and  Marsden  Ll 981b] ) . 
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§4.  An  Example:  the  simple  pendulum  coupled  to  two  oscillators 

We  illustrate  the  theory  developed  above  with  a  generalization  of 
our  earlier  two  degree  of  freedom  pendulum-oscillator  model  (Holmes  and 
Marsden  [1981b]).  Consider  a  simple  pendulum  linearly  coupled  to  two 
nonlinear  oscillators.  For  simplicity  we  assume  that  the  oscillators  are 
identical  (this  is  not  important)  and  that  their  Hamiltonians  can  be  ex¬ 
pressed  as  G(fc2  +  y^/2)  or,  equivalently,  in  action  angle  coordinates 
as 


G(I-),  i  =  1,2  , 

(4.1) 

with  Q(I - )  = 

^7  (1^*1*  0  for  Ii  >  0  , 

(4.2a) 

and 

n‘  (i • )  =  — ( i • )  t  o  , 

1  31?  1 

(4.2b) 

(cf.  eqns.  (2.3)  and  condition  (H2)  of  Section  2).  Elimination  of  either  1^, 

or  I2  by  reduction  is  then  possible.  For  definiteness,  we  shall  assume 

2  2 

that  I2  is  removed.  Our  assumption  of  the  form  G((x^  +  y..)/2)  is  merely 

for  computational  convenience,  since  more  ’realistic’  anharmonic  oscilla¬ 
tors  lead  to  Hamiltonians  G(T)  expressed  in  terms  of,  for  example,  elliptic 
functions  (cf.  Greenspan  and  Holmes  [1981]). 

The  system  to  be  studied  has  the  Hamiltonian 

He  »  p2/ 2  -  cos  q  +  G 1  Cl-j )  +  G(I2)  +  |  [(/HJ  sin  0j  -  q)2  +  (/2lJ  sine2  -q)2]. 

(4.3) 

The  unperturbed  orbits  in  the  homoclinic  manifold  are  given  by 

(q,p,91,62,I1,I2)  *  (±2  arctan  (sinh  t),  ±2  sech  t,  n(£ ^ ) t  +  6°  , 

nu2)t  +  e2,  a1  ,  iz)  , 


(4.4) 


i 
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where  h  >  1  is  the  total  energy  and  Tf  =  1  the  energy  of  the  homoclinic 
orbit  and  =  h^,  )  =  h  -  1  -  h^ .  Assumptions  H(l)  and  H(2) 

of  Theorem  3.1  are  therefore  satisfied  and,  in  view  of  (4.2a,b)  we  can  pick 


h  and  h^  so  that  the  nonresonance  conditions  necessary  for  application 


of  the  KAM  theorem  are  met.  To  check  the  final  assumption  we  compute  the 


Poisson  brackets  and  {F,H^}.  From  (4.3)  we  have 


{ 1 1  ,H^ }  «-!§-.»  -(^iy  sin  e1  -  q)  S2l]  cos  9] 


rc  U1 ,  3F  3H'  3F  3H'  .  n  r  ,  „  . 

{F’H  }  =  “apT "  3^“sF=  sinq*°  '  p[-{/2Ii  Sin  91  ‘  V 


-  (vflT  sin  0,  -  q)] 


=  p[^T7  sin  9-|  +  /2I2  sin  ©2  -  2q] 


Thus,  using  (4.4)  we  have 


Mi (0°,02,h,hi )  =  -[/2t^  sin(Q(£^)t  +  0°)  t  2  arc  tan(sinh  t)]< 


/2S^  cos(iT(£1  )t  +  0°)  dt  , 


^2^01 ,02,fl’^l  ^  *  -  2  sech  sin(n(£-|)t  +  0°)  + 


(4.5) 


si n(f2( Jl2)t  +  02)  +  4  arctan(sinh  t)]  dt.  (4.6) 


» 


Noting  that  the  integrals  of  products  of  odd  and  even  functions  vanish 
over  the  infinite  domain  and  taking  the  positive  branch  of  the  homoclinic 
manifold,  these  two  functions  become 


00 


Mi  =  2v2  ^ 

M2  = 


arctan(sinh  t)  sin(Q(jij  )t)dt*sin 

mm  OO 

00 

sech  t  cos(ft(£,)t  dt*sin  e° 

-CO 


+  2v2£2 


CO 

' 

sech  t  cos(ftU2)t)  dt*sin  e? 

^-00 


(4.7) 


For  brevity  we  write 

Q(  ^  )  =  u>j ,  q(  &2 )  =  ^2  ‘  '.4.6) 

To  evaluate  the  first  conditionally  convergent  integral  we  choose,  for 
computational  convenience,  the  limits  as  follows: 


lim 

N 


Ntt/co, 


-Ntt/uj-j 


arctan(sinh  t)  sin  i^t  dt  * 


lim  —  arctan(sinh  t)cos  w,t 
N-»|  “l  1 


Nit/ gj-.  /-Ntt/ ui.  me  w,  t 


r 

+  -L 


1  COS  1 


dt 


-Nir/oj,  1  -Nit^1  1  +  sinh  t 


■  if 

La 


sechct  cos  u^t  dt 


=  —  (itw.  )  cosech 

oj-j  1 


TTU)^ 

~~r 


(4.9) 
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The  final  integral  is  obtained  by  the  method  of  residues.  Similarly,  we 
have 


sech  t  cos  u).t  dt  =  tto)<  sech 
J  J 


mo 

2 


;  j  =  1,2, 


(4.10) 


and  thus 

Mi  =  2^/2^  cosech 


''mo-j  1 


sin  0^  , 


M2  =  2tt  co-j  sech 


TTlO-i 


sin  0“  +  /2jl2  ^2  sech 


1T00- 


sin  0°J. 


(4.11) 


One  obtains  a  similar  result  (with  an  appropriate  change  of  sign)  on  the 
negative  branch  of  the  manifold. 

We  therefore  find  zeroes  when  0°  =  kir,  0°  =  ju  for  all  integers  l,  k 
and  it  is  easy  to  check  that 


det  DM  = 


3M 

301 


aM 

301 


3M1  3M2 
302  36° 


=  ±871^2  /fAJl 2  cosech 

Trail 

2  J 

sech 

mz 

2 

.0  . .  _ 

t  0. 


-Ml 

@2  = 


(4.12) 


Thus  the  final  assumption  is  satisfied  for  suitable  choices  of  h  and  h^ 
and  we  have 


4.1  Theorem.  For  e  sufficiently  small  the  Hamiltonian  system  (4.3)  has 
a  set  of  2  dimensional  invariant  tori  of  positive  measure  each  of  whose  un¬ 
stable  manifolds  intersects  its  stable  manifold  transversely.  Moreover,  a 

finite  transition  chain  of  such  tori  T1,  ...,  t"1  can  be  chosen  such  that 

e  e 

Wu(Tj)  ^WS(j+1)  and  WU(TJ+1 )  TR  Ws(Tj);  j  =  1 . m-1 .  Thus,  orbits 

£  £  £  £ 
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can  be  found  which  pass  from  the  neighborhood  of  any  torus  ‘P  to  the 
neighborhood  of  any  other  torus  T£  in  the  chain.  This  situation  obtains 
on  every  energy  level  He  =  h  >  1. 

Remarks .  1.  Arnold's  [1964]  example  is  similar  to  ours  in  some  respects, 

but  he  employed  explicit  external  forcing,  taking  a  t-periodic  two  degree 
of  freedom  system  He(q,p,0,  I , t)  rather  than  a  three  degree  of ■ freedom 
autonomous  system.  This  perturbation  was  further  chosen  to  vanish  on  the 
tori,  so  that  the  perturbed  tori  lie  in  the  same  positions  as  the  unperturbed 
tori.  As  we  remarked  in  Section  3,  this  is  not  necessary  since  the  bracket 
<d(X  q,X  i )  =  {L0,!.1}  vanishes  on  the  unperturbed  isotropic  tori,  and  thus 

the  integral  of  the  Poisson  bracket  along  the  unperturbed  orbits  still  pro¬ 
vides  a  good  measure  of  the  separation  of  the  perturbed  manifolds. 

2.  Although  the  theorem  asserts  that  diffusion  occurs  on  every 
energy  level  h  >  1,  the  latitude  available  for  choice  of  h1  and  hence  for 
satisfaction  of  the  nonresonance  conditions  increases  with  h.  Thus  the  'suf¬ 
ficiently  small  c'  approaches  zero  as  h  -*■  1 . 
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§5.  Conclusions 

This  paper  and  its  companions  (Holmes  and  Marsden  [1981b],  [1981c]), 
address  the  general  question  of  perturbations  of  integrable  multidimensional 
Hamiltonian  systems.  A  particular  area  of  interest  is  the  development  of 
a  method  for  investigating  the  integrabil ity  of  the  perturbed  problem,  and 
for  providing  a  qualitative  description  of  orbits  in  phase  space. 

In  the  present  paper  we  have  combined  a  reduction  technique  with  a  vec¬ 
torial  version  of  Melnikov's  [1963]  method  to  establish  the  existence  of 
Arnold  diffusion  in  Hamiltonian  systems  with  at  least  three  degrees  of  free¬ 
dom.  This  in  turn  implies  that  the  system  is  non  integrable  in  the  classical 
sense:  there  are  no  analytic  integrals  other  than  the  total  energy.  The 
method  is  applied  to  the  specific  case  of  a  pendulum  coupled  to  two  non¬ 
linear  oscillators.  It  is  shown  that  the  stable  and  unstable  manifolds  of 
non-resonant  tori  that  survive  under  a  small  perturbation  intersect  trans¬ 
versely.  We  briefly  discuss  how  this  enables  points  in  phase  space  to  dif¬ 
fuse. 
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